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Conformal Mapping on Rough Boundaries 
II: Applications to bi-harmonic problems 
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We use a conformal mapping method introduced in a companion paper Q to study the properties 
of bi-harmonic fields in the vicinity of rough boundaries. We focus our analysis on two different 
situations where such bi-harmonic problems are encountered: a Stokes flow near a rough wall and 
the stress distribution on the rough interface of a material in uni-axial tension. We perform a 
complete numerical solution of these two-dimensional problems for any univalued rough surfaces. 
We present results for sinusoidal and self-affine surface whose slope can locally reach 2.5. Beyond 
the numerical solution we present perturbative solutions of these problems. We show in particular 
that at first order in roughness amplitude, the surface stress of a material in uni-axial tension can 
be directly obtained from the Hilbert transform of the local slope. In case of self-affine surfaces, 
we show that the stress distribution presents, for large stresses, a power law tail whose exponent 
continuously depends on the roughness amplitude. 

PACS numbers: 02.70.-c, 46.30.Cn, 47.15.Gf 
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I. INTRODUCTION 

In a companion paper we have presented a confor- 
mal mapping technique that allows us to map any 2D 
medium bounded by a rough boundary onto a half-plane. 
This method is based on the iterative use of FFT trans- 
forms and is extremely fast and efficient provided that the 
local slope of the interface remains lower than one. When 
the maximum slope exceeds one this algorithm, similar in 
spirit to a direct iteration technique well suited to circu- 
lar geometries 0,13, can no longer be used in its original 
form. Underrelaxation'^ suffices however to make it con- 
vergent for boundaries having large slopes. Beyond the 
determination of a conformal mapping for a given rough 
interface we have also shown in Ref. how to gener- 
ate directly mappings onto self-affine rough interfaces of 
chosen roughness exponent. Self-affine formalism is an 
anisotropic scaling invariance known to give a good de- 
scription of real surfaces such as fracture surfaces 0,0, 13 ■ 
This statistical property of fracture surfaces is of great 
interest in the study of friction or transport processes in 
geological faults 0, 13- 

Building a conformal mapping is entirely equivalent to 
solving a harmonic problem with a uniform potential (or 
field) condition on the boundary. We used extensively 
this property in Ref. 1] to study stationary heat flows in 
the vicinity of a rough boundary and we focussed on the 
case of self-affine surfaces where we were able to com- 
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pute the exact correlation between local surface field and 
height profile. We gave also a special emphasis to the 
problem of the location of the plane interface equivalent 
to the rough one at infinity. It turned out that the con- 
formal mapping technique provides a very direct means 
of computing the shift between the plane equivalent in- 
terface and the mean plane of the rough interface. 

The range of applications of this first study naturally 
covers fields where Laplace equation appears: electro- 
statics, concentration diffusion, antiplane elasticity... In 
this second paper, we are specifically concerned with the 
case of bi-harmonic problems. The method we propose 
leads to the solution of the bi-harmonic field through the 
inversion of a linear system, as most other alternative 
numerical approaches (e.g. boundary elements method, 
spectral method...) but in contrast to the latters, the lin- 
ear system to invert is naturally well conditionned and of 
rather modest size (N equations for N Fourier modes in 
the conformal mapping method) in contrast with direct 
spectral methods (4N unknowns). 

Moreover, following the first step of our algorithm ana- 
lytically allows to obtain systematic perturbation expan- 
sion results. 

After recalling our main results about conformal map- 
ping in the first section, we deal successively with two im- 
portant examples of bi-harmonic problems: Stokes flows 
and plane elasticity. The second section is thus devoted 
to the study of a stationary Stokes {i.e. low Reynolds 
number) flow close to a rough boundary. We shall also de- 
velop in this section the paradigm of the equivalent "no- 
slip" plane interface. In the third section we point out the 
problem of the stress distribution in a two-dimensional 
material bounded by a rough boundary and submitted 
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FIG. 1: A schematic illustration of the mapping fl which maps 
the semi-infinite plane 2? onto the domain limited by a rough 
interface £. 

to a uni-axial tension. The study of such situations is 
of particular interest for computing the influence of the 
roughness on the rupture probability law of brittle mate- 
rials. We show that this problem is formally identical to 
that previously solved for the Stokes flow. We pay fur- 
ther a special attention to the case of self-affine surfaces 
and we present in this section numerical results that sug- 
gest that for large stresses, the surface stress statistical 
distribution law presents a power law behavior. 



our study to mappings Q that are bijective holomorphic 
functions, i.e. only depends on the variable z and not 
on its conjugate z — x — iy. The transformations associ- 
ated with such functions are said to be conformal in the 
sense that they preserve locally the angles. Let us now 
take advantage of the semi-infinite geometry we have to 
deal with. The two domains we consider are very similar 
apart from the region close to the boundary. Far from 
this one, Q is essentially the Identity and we can write: 

n{z) ^ z + uj{z) , (1) 

where the perturbation lu decreases with depth and takes 
non negligible values only in the close vicinity of the in- 
terface. In the following, we consider periodic interfaces 
in order to minimize edge-effects. Without loss of gener- 
ality, let us choose 27r to be the lateral period. A natural 
form of Q is then: 

oo 

n{z)^z + u}{z) = z + ^ujke-'^''''' . (2) 

k=Q 

where uj is expanded on a basis of evanescent modes. 



II. CONFORMAL MAPPING ON ROUGH 
BOUNDARIES 

We recall here the essential results described in Ref. 
0. The aim of this section is to conformally map a 
half-plane onto a two-dimensional domain bounded by 
a rough interface. We first recall briefly how to build 
a conformal mapping well-suited to a given rough inter- 
face. This problem can be written under form similar as 
that of a "Theodorsen problem" Q ; in the semi- infinite 
geometry we deal with, it can be solved with an iterative 
algorithm using fast Fourier transforms 0. The second 
part of this section is focussed on the specific case of self- 
affine interfaces. It turns out, indeed, that very simple 
constraints on a conformal mapping allow to generate di- 
rectly a two-dimensional domain bounded by a self-affine 
interface of chosen roughness exponent. This property is 
of particular interest in statistical studies. We could thus 
establish in Ref.0 the correlation between the norm of 
a harmonic field at a self-affine interface and the height 
of this latter. 



A. Notations 

As illustrated in Fig. ^ we are seeking for a map- 
ping from a half-plane onto a two-dimensional domain 
bounded by a rough interface. In the following we place 
our study in the framework of the complex analysis. We 
consider then the lower half-plane T) whose complex co- 
ordinate z = X + iy is such that Smz < and the two- 
dimensional domain £ bounded by the rough interface 
d£, we call w = u + iv the complex coordinate in £. We 
are seeking for a mapping Q from V onto £. We restrict 



Computing the mapping for an imposed 
interface 



We consider a single valued interface d£. Let h be the 
real function giving the interface geometry, for all point 

w = u + iv of d£, 

V = h{u) . (3) 
The mapping fimction Q is such that Q {dV) — d£, i.e 



ikx 



\k=0 
oo 



(4) 



where it = x + 5Re j 



-ikx 



\k=0 



The first equation is here very close to a Fourier trans- 
form except that we have h{u) instead of just h{x) in the 
first term. This formal proximity can be used to build an 
iterative algorithm. For sufficiently small roughness, we 
can see from the second equation that x is an approxima- 
tion of u at zeroth order. A direct Fourier transform of 
the profile h{x) allows then a first approximation {lo^^} 
of the coefiicients usk- The latter can be used to correct 
the previous approximation of u{x), using Qj) gives then 
a following estimation of the non-uniform sampling u{x) 
and of the coefficients ujk via the Fourier transform of 
h[u{x)]. In appendix A, we give a more formal presen- 
tation of this algorithm. It turns out that this iterative 
technique converges provided that the maximum slope of 
the profile remains below unity. The technique can be 
used for any rough single-valued interface. For profiles 
whose maximum slope exceeds one, the algorithm can be 
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made convergent with slight modifications such that the 
use of under-relaxation techniques. We refer the reader 
to Ref. Pi for more details on the convergence and the 
stability analysis in this specific framework. Extensive 
studies of this techiiique in case of cicular geometry are 
available in Ref. 0, [J. 



C. Conformal mapping on self-afflne interfaces 

As pointed out above, the algorithm we have just de- 
scribed is suited to any rough interface. It is in particu- 
lar possible to build conformal mappings associated with 
self-affine interfaces. The latter are defined by their scal- 
ing invariance properties: an interface described by the 
equation y = h{x) is said to be self-afiine if it remains 
(statistically) invariant under the transformations 



Ax , 
AS 



(5) 



for all values of the real parameter A. The exponent C is 
called the "Hurst" or roughness exponent. It is charac- 
teristic of the scaling invariance. From this property, we 
derive easily that 



{{h{x)-h{x + 5)f) ^C^5^^ 



(6) 



where C is a prefactor. A simple Fourier transform gives 
then the power density spectrum of the rough self-affine 
profile. 



P[k) cx k: 



-1-2C 



(7) 



When using the algorithm previously described to map a 
self-affine interface, the first guess for the mapping func- 
tion coefficients ojk is thus such that: 



ujk = 2iak oc k' 



Xk 



(8) 



where the ak are the coefficients of the Fourier transform 
of the profile and Xk are k independent random variables. 
It turns out that this power law behaviour is not altered 
by the following steps of the algorithm. In a symmetric 
way, we can impose, without any further restriction, the 
LUk to follow a power law and have a look at the interface 
generated. We choose thus 



(9) 



where e/c are random Gaussian variable with mean and 
unit variance for the real and imaginary part indepen- 
dently. We have however to note that nothing prescribes 
a priori that the function obtained is bijective. From the 
parametrical expression of the interface d£, we can write 



5Re 



dx 



{x + iO) 



= 1 



E 

. k 



and to guarantee that d£ remains single- valued we have 
to choose amplitudes lower than the threshold value 
Amax where 



A — 



-1 



(11) 



(10) 



We checked numerically (see Ref. Ij that the power spec- 
trum of such synthetic profiles was indeed a power law 
with the expected exponent. 



III. STATIONARY STOKES FLOW IN THE 
VICINITY OF A ROUGH WALL 

The Stokes equation describes fiuid flows at low 
Reynolds numbers. We address in this section the prob- 
lem of a stationary Stokes flow in the vicinity of a rough 
boundary. The study of such flows can be of great tech- 
nological interest in the case of convective transport pro- 
cesses one can think of problems of surface depo- 
sition or erosion. In the same spirit, the occurrence of 
recirculating eddies can render very difficult the decon- 
tamination of a polluted surface ; contaminants parti- 
cles can be captured by diffusion in a cavity and remain 
trapped in it for an arbitrary long time. From a more 
fundamantal point of view, in experiments consisting of 
tracing particles passively advected in a flow, the same 
process can lead to non gaussian statistics of the arrival 
time of tracer particles. 

We consider a semi-infinite 2D geometry with periodic 
lateral boundary conditions and a unit shear rate at in- 
finity. Far from being specific, the results obtained in 
this context can easily be extended for any case of sta- 
tionary shear flow in the framework of a double scale 
analysis. Considering a shear flow of shear rate 7 and an 
interface of typical boundary e, where the the velocity 
and pressure fields U and p obey the usual Navier Stokes 
equation, we can, following Richardson 11], define an 
inner problem where the reduced non dimensional vari- 
ables obey, at first order in e, a simple Stokes equation 
in a semi-infinite geometry. In the following we present 
a solution of this Stokes equation that can be rewritten 
as a biharmonic equation for the stream function. This 
solution only requires the knowledge of a conformal map- 
ping n from the lower half plane V onto the actual space 
£ bounded by the rough interface d£ and the inversion 
of a well conditioned linear system ; it can thus be ap- 
plied to any single-valued rough interface. We focus this 
brief study on the problem of the determination of the 
location of a plane boundary equivalent to the rough one 
at infinity. This problem is equivalent to the one of the 
replacement of the no-slip condition on a rough interface 
by a back flow condition (to be determined) on the mean 
plane. We compare our results with those of Tuck and 
Kouzoubov fl3| who developped in the actual space a 
method similar to ours in spirit. Recent results about 
Stokes flows near rough boundaries can also be found in 
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references ^2 ^3 ; in most of them the Stokes equa- 
tion is solved using boundary elements methods (see for 
instance the review of Pozrikidis [THjl). 



A. General solution 

We adress here the problem of a unit shear Stokes flow 
in the vicinity of a rough boundary. Let us call ^'(it;) the 
stream-function associated to the veolociy field b in the 
actual space E. We have by definition 



dv 



du 



(12) 



In a stream function formalism, the Stokes equation is 
reduced to a simple bi-harmonic equation. Taking into 
account the boundary conditions i.e. no slip on the in- 
terface and unit shear rate at infinity, '^{w) has to be 
solution of the following problem: 



Vi^{w) =0 in 
Vu,*(w) = on 
*(w) ^v"^ as 



E 



(13) 



The essential difficulty obviously lies in the no slip con- 
dition on the interface ; the use of a conformal mapping 
allows us to build an equivalent problem with a much 
easier boundary condition, the new interface being plane 
instead of rough. Let us associate to the stream-function 
^ in the actual space £ the real potential <& in the half- 
plane V: 



$(z) = ^ oVL-^{w) 



(14) 



where J7 maps V onto £. We can thus define the following 
equivalent problem in the new geometry: 



Vt*(i«) 
V„^'(u;) 
*(w) 

.WW 

V^$(z) = 





on d£ 



as V 



= 



(15) 



<i>(z) 



on dV 

as y ^ 



In the case of a simple harmonic equation, building the 
conformal mapping gives immediatly the complete so- 
lution; this is unfortunately no longer true in the case 
of a bi-harmonic equation. One can see that the origi- 
nal equation is changed into a linear equation with non- 
constant coefficients. The latter equation is directly re- 
lated to the mapping function fl. We show in the fol- 
lowing that this difficulty can be circumverted. Let us 
recall that in addition to the above described conditions, 
the new potential $ has to be real and 27r-periodic in 
X. Besides, the boundary condition can be made simpler 
taking into account that the interface is now plane. <I> 



being defined apart from an additive constant, we can 
write that it obeys: 

f $ =0 ondV 
V.$ = ondV ^ a$ ^0 ^^^Q^ (16) 

The bi-harmonic potential $ can always be written in 
terms of two holomorphic functions F and H such that: 

$(z) = n{z)F{z) + F{z)n{z) + H{z) + i?(z) . (17) 

In the following we split both functions into a purely pe- 
riodic part (denoted by the index p) and a non-periodic 
part. Taking into account the desired behaviors at infin- 
ity provides: 
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F{z) ^ -z + Fp{z) , 
H{z) =. -L^ + Hi{z) + Hp{z) 



with 



n>0 



(18) 



(19) 



and Hi is z times a x-periodic function. The lateral peri- 
odicity of $ forbids the occurrence of terms proportional 
to polynoms of the real variable x, hence: 



Hi{z)^-z[Fpiz) + ^^iz) 



(20) 



and $ becomes then: 
1 



n n 

+ 2iyY Un + l^n 

n ^ 



n>l p 



n>l p 



+ ^(l.„/„ + cJ„/„) e2»y . 

n 

(21) 

We can deduce from the latter expression the partial 
derivative of $ with y: 

H (^) = y + E„ nKe-"^^ + E„ "7^«e™" 

+2iy « (/« + l^n) e"'~ - 2iy E„ n (/„ + Jl^„) e" 



- E„>1 Ep + 2p) {fpUJn+p + ^pln+p) 

- E„>1 Ep + 2p) {uJpfn+p + Ip^n+p) 

-E„2n U7„ + cJ„/„) e2»^ . 



^inx ^(n+2p)y 
~inx (n+2p)y 



(22) 
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The holomorphic functions fip, Fp and Hp being 2-k- 
periodical, they can be developed using the basis of the 
functions {e~*'^^}. Besides, the interface dV being simply 
the X-axis, the restriction of these holomorphic functions 
to the boundary can be written in {e^*'"'^}. The bound- 
ary condition problem can then be written by cancelling 
out the projections of $ and dy^ on the functions- vectors 
|g-ifca;j._ Keeping only the first N modes, we obtain 2N 
equations which allow to obtain the 2N components /„ 
and /i„. The projection of the boundary condition on the 
function- vector 1 (fc = 0) gives first: 



■ X! (^»/« + ^"/") + 2i ( /o + 



2i /o + -too = 



(23) 

We have then for the function- vectors {e *'^^} with fc ^ 0: 

' $ : ^ {uJpfk+p + Ipi^k+p) +hk^O , 
p 

■ + 2p) ippfk+p + Ip^k+p) + khk (24) 



2* [Ik + g^fc j - . 

We can thus write the following linear system: 
1 



fk + gto/c +i^P i^pfk+p + fp^k+p) = 
p 

^fc = - X! i^pfk+P ~ fp'^k+p) ■ 



(25) 



The coefficients {/&} are solutions of the first equation, 
a. N X N linear system. The coefficients {hk} are easily 
deduced from the {fk}- Once the conformal mapping is 
known, the numerical resolution of the Stokes equation 
just requires the inversion of the linear system. Let us 
note that the latter system is easy to invert, which would 
not have been the case if we have written the bound- 
ary condition in the ui-space. The latter method was 
recently used by Tuck et Kouzoubov 0. In the cited 
reference, they use expansions in a basis of {e~'^^cos(A;a;)} 
and write a linear system discretizing the rough interface 
in N points. If this method is efficient in the small slope 
limit, it becomes however unpracticable for large slopes. 
The procedure requires thus the numerical inversion of a 
matrix consisting of terms of order e^^"^ where A is the 
roughness amplitude and this becomes practically diffi- 
cult or unprecise as soon as the product NA increases. 
The conformal mapping avoids such numerical difficul- 
ties since the boundary is plane in the equivalent domain. 
One can see on figures |21 and 13 maps of the stream func- 
tion for a stationary Stokes flow for two rough surfaces 
identical up to a dilation of factor 4 in the vertical direc- 
tion. The stream lines closely follow the smooth interface 
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FIG. 2: Stream lines of a Stokes flow along a rough surface. 



-1.0 - 




FIG. 3: Stream lines of a Stokes flow along the same rough 
surface as above but four times rougher. We observe that a 
recirculation flow appears in the deepest cavity. 



(Fig. 121) while an eddy appears in the largest depression 
of the roughest interface (Fig. 



B. Equivalent plane boundary 

We now introduce the notion of an equivalent plane 
no-slip boundary in the framework of a stationary Stokes 
flow. Our aim is here to replace the no-slip condition on 
the rough boundary by a no slip condition on an equiva- 
lent plane boundary, the stream function remaining un- 
changed at infinity. Let us recall that like in the case of 
a rough electrode which was studied in Ref. [J, noth- 
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ing prescribes the plane equivalent boundary to lie at 
the average height of the rough one. In harmonic prob- 
lems, the dissymetry between peaks and cavities due to 
point effects leads the equivalent plane boundary to be 
shifted towards the peaks. We shall see in the following 
that the same conclusion holds in the case of the sta- 
tionary Stokes flow. An illustration of this disymmetry 
naturally emerges with the occurrence of little eddies in 
pronounced depressions of the rough boundary. We have 
chosen in this text to develop the paradigm of no-slip 
plane equivalent boundary but it is naturally also pos- 
sible to consider a plane boundary fixed at the average 
plane with a roughness dependent slip boundary. The 
conclusion previously mentioned about the place of the 
equivalent boundary (nearer from the peaks than from 
the cavities) is then expressed by a reversal flow condi- 
tion. We show in the following that the conformal map- 
ping method gives a natural way to compute the vertical 
shift of the plane equivalent boundary. We compare these 
results with a perturbative solution that can be directly 
computed from the Fourier coefficients of the interface. 



1. Conformal mapping approach 

As already mentioned, the stream function ^'(u,?;) of 
the Stokes flow is entirely determined by the following 
conditions: 



Vt*(w) =0 in 
V„*(w) = on 



as V — oo 



(26) 



If we now return to the solution obtained via conformal 
mapping, we have: 



J2n>l J2p {fp^n+p + 



Jnx {n+2p)y 



+ E„>i Ep K/n+P + 7p^n+p) e-™-e("+2f)J' 
+ E„ K7„ + cJ„/„)e'"^ • 

(27) 

By construction, ^{z) is such that ^'(ly) is bi-harmonic 
in £ and fulfills the no-slip boundary condition at the 
rough interface. Let us now build the stream function 
^eq associated with a plane interface located a,t v — H; 
we have immediatly: 

^eg{w)^^{v~H)\ (28) 

and its associated function in the half-plane V is 

<i>,g{z) = l['^mn{z)-Hf , (29) 



that becomes at infinity: 

$(z) = ^ 1 y cj^ ^ 8j^) ^ as y -oo , 

^eq{z) = \y^ + y i'^M^o) -H] + 0(1) as y ^ -oo . 

(30) 

Taking into account Ea. (|23|l which specifies the expres- 
sion of /o, the identity between $ and <f>eij defines the 
value of H 

H = 5to(cjo) + 2 ^ n (a;„/„ + ujnf„) . (31) 



2. Perturbative approach 

Let us consider an interface of amplitude, say A, with 
characteristic length A such that the profile is statistically 
symmetrical. When we seek for the location of the plane 
equivalent electrode, we expect two different behaviors 
i) a dependence on A^/X in case of small amplitude or 
low spatial frequency ii) a linear dependence on A in case 
of large amplitude or high spatial frequency. The latter 
behavior comes directly from the fact that the equiva- 
lent plane reaches at maximum the level of the highest 
peaks. In case of small slopes, it is easy to show that the 
correction from the average plane is of order A"^ /X. The 
deviation H is naturally normalized by the amplitude of 
roughness A. The ratio H/A has then to be a function of 
the two only characteristic lengths of the system, A and 
A and can be expanded in the limit of small slopes: 

H ,,A, A /aV ,^fA\^ , , 

^=0(^)^ao + «r^ + a.y . (32) 

A simple symmetry about the mean plane has to leave H 
unchanged, since this symmetry is equivalent to a trans- 
formation of A into —A, oq and a-i have to be zero and: 



R = ar 



4! 

A 



A^ 
A3 



(33) 



A detailed perturbative analysis can be built in the 
case of a simple sine interface 10]. Writing a perturbative 
solution in the conformal mapping formalism allows us to 
deal with any rough interface. Following (|21|l we have: 



$(z)= ^y^ + Ay^m (Ppiz) + ^cu{z) 



+2'!fte [cjiz)Fp{z) + Hp{z) 



(34) 



By construction, cj, Fp and Hp are of order A, A being the 
roughness amplitude. At first order the no-slip boundary 
condition becomes: 



$(x) = 
9$ 



dy 



(^H^p^\x)'^ = 0, 
4am(^F(i)(x) + lc.(i)(a;)) (35) 

. 



+25Re I ^^{x) 



dy 
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The holomorphic functions being bounded, the resolution 
of these Hilbert problems gives immediatly: 



i/W(z)^0 



p 



(^) 



1 



(36) 



Using the iterative algorithm briefly presented in section 
II (see Appendix for details), a first order approximation 
of the coefficients LOk is 



.hk_ 

'' N 



(37) 



where {h} is the 2N discrete Fourier transform of the real 
function h associated with the rough interface. Using the 
expression of $Jm(cjo) derived in section V-B of Ref. 0, 
we can write: 



H 



1 

'iV2 



(38) 



fc>0 



This result is consistent with the one proposed in Ref. 
for a backflow slip condition on the mean plane of 
the interface. In the particular case of a pure sine profile 
of amplitude A and wavelength A, we recover: 



H 



-2tt- 



4! 

A 



(39) 



It has to be noted these first order results are exactly 
identical to those obtained in the case of a rough elec- 
trode (up to a factor 2) despite the fact we had here to 
solve a biharmonic equation instead of a simple harmonic 
one. On figures 01 and [3 we have plotted both results of 
the perturbative solution and of the conformal mapping 
calculation in case of a sine interface and of a self afRne 
interface of roughness exponent 0.8 built with 64 Fourier 
modes. We check that the perturbative calculations cor- 
rectly fit the results for small slopes (up to 0.5). For 
larger slopes, the perturbative expression overestimates 
the deviation whose behavior becomes progressively lin- 
ear. Note that our numerical method allowed us to reach 
local slope values up to 2.5. This maximum slope can be 
increased by using more Fourier modes in the mapping 
function (we used 256 modes in the present calculation). 
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FIG. 4: Shift from the mean plane of the plane equivalent 
boundary for a Stokes flow in case of sine interface of varying 
amplitude A. The dotted line corresponds to the second or- 
der perturbative expansion and the symbols to computations 
using conformal mapping. The dashed line has a slope -1. 
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FIG. 5: Shift from the mean plane of the plane equivalent 
boundary for a Stokes flow in case of self-afBne interface of 
varying amplitude. We use in abscissa the maximum local 
slope Smax- The dotted line correspond to the second order 
perturbative expansion and the symbols to computations us- 
ing conformal mapping. The dashed line has a slope -0.06. 
The surface has been built with 64 Fourier modes and we 
used 256 modes in the solution based on conformal mapping. 



IV. ELASTICITY 

In this section, we analyse the effect of a slight sur- 
face roughness on the stress distribution in an elastic 
medium. We emphasize here the case of a semi-infinite 
material submitted to a uni-axial tension. Although very 
elementary, this simple model illustrates an effect which 
has been suggested to be responsible for the mechanical 
strength of glass fibers. Recent experimental results 
suggest that the nanometric roughness at the surface of 
glass fibers of diameter a few micrometers could be re- 
sponsible for the decrease of the tensile resistance by a 



factor of about 5. In uniaxial tension, the resistance of a 
fiber is directly related to the distribution of maximum 
positive principal stresses. Using a simple perturbative 
expansion, we show that in the limit of small slopes, the 
surface stress can be directly computed from the Hilbert 
transform of the local slope. 

The WeibuU law 0| usually gives a correct description 
of failure statistics in a wide range of brittle materials. 
This phenomenological law is partly based on the iden- 
tification of the weakest link in the material (whose size 
is supposed to be that of the smallest defects). We will 
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see in the following that for a self-afRne surface (such as 
measured for glass fibers by AFM techniques 16] ) the 
statistical distribution of tensile stresses at the bound- 
ary displays a power law behavior, which naturally im- 
plies the validity of WeibuU statistical failure distribu- 
tion. Moreover the exponent of this power law, i.e. the 
WeibuU modulus, continuously depends on the roughness 
amplitude. 



A. General solution 

In plane stress or plane strain conditions, the stress 
tensor [a] can be completely represented by a unique real 
function, named the Airy function: 



^ uv 



where < 



dudv 

9^ 



(40) 



This form directly comes from the stress balance in ab- 
sence of external forces i.e. div[cr] = 0. In the framework 
of 2D elasticity in an isotropic medium, the Airy function 
obeys: 



AA* = in 5 



(41) 



The stress tensor being computed from two successive 
derivations of the bi- harmonic function the latter is 
only defined apart from a linear function in u and v. 
In the following, we consider free boundary conditions 
[(j]n = and we impose a uni- axial tension at infinity. 
Let n denote a unit vector normal to the interface, we 
have: 







1 




on dS , 

as y ^ - 



(42) 



The Airy function 5* is thus such that 



n. 



"9^ 



n 

" dudv 



= 



on 88 
on d£ 



(43) 



4' 



1 2 

as V 



At any point of the interface d£, it is possible to give a 
parametric representation of the tangential and normal 
vectors t and n: 



t{w) 



VL' (x) 

\n' {x)\ 



n {w) 



in' (x) 



(44) 



The boundary conditions at the interface can be rewrit- 
ten: 



[crjn = on d£ 

^ (t ■ V^) Vu,* =0 on d£ 
^ V^* = Cst. = on d£ . 



(45) 



We can choose the constant to be zero since the Airy 
function is only defined apart from an affine function. 
This leads to: 



ondD 

<^ = on dD 



(46) 



It turns out that the boundary condition at the interface 
is exactly the same as the one we have encountered for 
the no-slip condition in a Stokes flow. The Airy function 
is thus identical to the above derived stream function for 
Stokes flow. 



B. Surface stress distribution 

1. Perturbative approach 

The normal stress being zero at the interface, the first 
order expression derived in the previous section gives us 
the following result for the principal (tangential) stress 
at the interface: 



= l-25Re pi)'(x) 

k—7i J 'T"' 

^ \ky '2n 

= 1 - 2H[ip']{x) , 



(47) 



where Ti,[...] stands for the Hilbert transform oprerator on 
the real axis. At first order, we find that the tangential 
surface stress deviation from its mean value is propor- 
tional to the Hilbert transform of the local slope. Fig. |^ 
andClgive two examples of stress distribution on surfaces 
of repectively maximum slope 0.1 and 0.4. In both cases, 
we can see that the stress fluctuations are much greater 
than the height fluctuations (which have been dilated by 
a factor 5 in the flgures). Note that (taking into account 
the dilation of the height profllc) the stress fluctuations 
are very large compared with those of the height. In 
the context of rupture, the relevant parameter being the 
maximum stress, one can understand that a very modest 
roughness can be responsible for a dramatic decrease of 
the material resistance. 

We observe a good agreeement between the stress pro- 
flle computed by conformal mapping and the perturba- 
tive result (related to the Hilbert transform of the slope) 
for the smoothest interface but this is no longer the case 
for the roughest one especially in the area where the cur- 
vature is large. Higher order perturbative terms are thus 
necessary to recover the actual stress. 
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FIG. 6: Stress and height profiles on the rough interface of 
a 2D medium submitted to uni-axial tension. The bold line 
represents the height profile (of maximum slope 0.1) dilated 
by a factor 5. The solid line gives the stress distribution 
obtained by the complete conformal mapping computation, 
the dotted line is a first order expression of the stress which 
is directly obtained from the Hilbert transform of the local 
slope of the interface. 



-0.10 -0.05 0.00 0.05 0.10 

log,o(o) 

FIG. 8: Surface stress distribution for self-affine surfaces of 
roughness exponent = 0.8 and of roughness amplitudes 
A = O.Obemax, tmax is the amplitude such that the maxi- 
mum local slope is equal to 1. The self-affine interfaces have 
been built with 32 Fourier modes and the results averaged on 
1500 surfaces. We can check that the distribution (symbols) 
is well fitted by a parabola, which shows that the stress is 
log-normally distributed. 



1.0 




FIG. 7: Same figure as above with an interface four times 
rougher. One observes that in large curvature areas, the first 
order approximation does no longer suffice to represent pre- 
cisely the local stress. The height profile (of maximum slope 
0.4) dilated by a factor 5. 



2. Statistical results 

Let us now turn to the study of stress distributions on 
self-affine surfaces. The latter are designed to present a 
Gaussian height distribution i. e. their Fourier coefficients 
are hk — Aekk^'^~'' where ek is a Gaussian random vari- 
able of zero mean and unit standard deviation. In case 
of small slopes, the validity of the first order perturba- 
tive result suggests thus that the stress distribution is 
Gaussian. On Fig. |Slwe have plotted in log-log scale the 
surface stress distribution obtained for 1500 self-affine 



surfaces of roughness exponent C = 0.8 and of maximum 
slope 0.05, one can check that, as expected, the distribu- 
tion is well fitted by a parabola. 

In case of larger slopes we have seen that the agreement 
between the first order perturbative results and the com- 
plete computation becomes poorer. We can then wonder 
if the log-normal behavior of the distribution is preserved. 
On Fig. |5| we have plotted the surface stress distribu- 
tion for four self-affine surfaces of roughness exponent 
( = 0.8 and of respective maximum slopes 0.2, 0.4, 0.6 
and 0.8. These results were obtained by averaging the 
data obtained from 1000 different surfaces each defined 
with 64 Fourier modes. We can see a clear power law 
like behaviour for large stress amplitudes. The slopes we 
can measure are very dependent on the roughness ampli- 
tude. The interpretation of these new numerical results 
requires a perturbative analysis that we have not devel- 
oped yet for this bi-harmonic problem. We have however 
performed a similar analysis in the case of a harmonic 
field on self-affine interfaces It turned out that in a 
similar fashion as the one presented above, the field dis- 
tribution law present a power law tail with an exponent 
T oc A~'^£-^~'' where A is the roughness amplitude, £ the 
spatial lower cut-off of the self-affine domain and C the 
roughness exponent. Calling g the logarithm of the field, 
the latter result was derived showing that the reduced 
variable 

fAig) - ( + ^Kg - l) /KA (48) 

follows a Gaussian distribution {K ~ 2 for harmonic 
problems). Using K — 0.25, we can check indeed (see 
Fig. I10|l that all data obtained from our calculations 
collapse on a same parabola in log-log scale. These nu- 
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FIG. 9: Surface stress distribution for self-afBne surfaces of 
roughness exponent C, = 0.8 and of roughness ampUtudes 

A = 0.2€max (A), OAe^ax (o) , OSe^ax (H) and O.Stmax (o), 

^max is the amphtude such that the maximum local slope is 
equal to 1. The self-afRne interfaces have been built with 64 
Fourier modes and the results averaged on 1000 surfaces. For 
each distribution, a bold line shows the power law behavior 
obtained for large stresses. 




-3.0 I ' ' ' ' ' ' ' ' ' ' '■ 1 

-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

fA[log,„(a)l 

FIG. 10: Same distributions as in the previous figure in 
the reduced variable fA{g) = (\/l + '^Kg — l) /KA where 
g = logj^Q(a) and K — 0.25. The data collapse onto a simple 
parabola, which shows the Gaussian character of the distri- 
bution of /a ((?) 

merical results indicate that a same scaling applies for 
both harmonic and bi-harmonic fields on self-affine sur- 
faces. A detailed analysis of these statistical properties 
will be addressed elsewhere. 



V. CONCLUSION 

After introducing a conformal mapping technique that 
allowed for a detailed study of harmonic field in the vicin- 
ity of rough boundaries [l| , we have extended in this pa- 



per the use of this method to the study of bi-harmonic 
fields. We have given a general solution to problems such 
that the Stokes flow over a rough surface and the stress 
distribution in a medium (bounded by a rough inter- 
face) in uni-axial tension. Besides the knowledge of the 
mapping function (obtained using a simple iterative al- 
gorithm), this solution only requires the linear inversion 
of a well conditioned matrix. The determination of the 
mapping function being only limited by the maximum 
value of the local slope at the interface, the method is 
well suited to any kind of single-valued interface. As an 
illustration, we have thus presented results of Stokes flow 
over self-afSne boundaries whose maximum slope reach 
2.5. In the same context of a Stokes flow over a rough 
boundary, we pay a special attention to the determina- 
tion of the location of an equivalent no-slip plane inter- 
face. The conformal mapping method gives way to a very 
natural determination of this quantity. A simple pertur- 
bative solution allowed us to retrieve for it a result pro- 
posed by Tuck and Kouzoubov ^3 • t;he context of the 
plane elasticity, the same perturbative result has allowed 
us to show that, in the limit of small slopes, the sur- 
face stress distribution was directly related to the Hilbert 
transform of the slope of the interface. This very simple 
result could be used e.g. for the evolution of a stress- 
corrosion front. The analysis of statistical results for the 
principal stress on self-afRne surfaces has shown moreover 
that the large stress distribution presents a power law tail 
whose exponent continuously depends on the roughness 
amplitude. Such results could be applied in the context 
of glass fiber rupture statistics to provide a fundamental 
basis for the WeibuU law that is known to describe ac- 
curately the rupture statistics. A realistic description of 
these stress distributions require a second order pertur- 
bative analysis which will be presented in a futher study. 
We have, however, recent ly p roposed such an approach 
in case of harmonic fields 18] where distributions of the 
same kind have also been found and justified in two and 
three dimensions. 
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Appendix: Iterative algorithm 

In this appendix we present the iterative scheme al- 
lowing us to build the conformal mapping Q{z) — z + 
uJkG^^^^ for any given single valued rough interface. 
As exposed in II. B, we are seeeking for a mapping func- 
tion such that 17 (dV) — d£, where d£ is the rough 
interface and dV is the x-axis. The single valued in- 
terface d£ being associated to the real function h, this 
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condition is rewritten: 



UJkC 



-2i7:kz 



(49) 



where it = x + Re ( 



LUke 



-2inkz 



We build the algorithm in the following way: the in- 
termediate quantities appearing at the fcth iteration are 
labelled with a superscript (fc), all functions are decom- 
posed over a set of 2n discrete values, the number of 
Fourier modes is thus limited to 2n. We first introduce a 



(k) 

series of sampling points u) with j — 0, 



(0) 



is initially set to an arithmetic series u 
sampling of h(u) by the u^*''' gives the array: 



71—1 which 
= jn. The 

(50) 



The discrete Fourier transform of this array is the com- 
plex valued array 



(51) 



for — n < j < n. The latter is shortly written as: 

a(^) = ^[/jW] , (52) 

where denotes the Fourier transform, which will be 
chosen as the Fast Fourier Transform (FFT) algorithm, 
thus imposing that n is an integer power of 2. The inter- 
mediate mapping w'^'^^ is computed from the a*^*^^ as: 



,(fc) _ 

u„ - 
,(fc) _ 



{^/n)af^ 

(i/2n)4''^ 





for j > 



for j < 



(53) 



The latter form is obtained from the identification of 
Ea. (|49b ) and the definition of a^''\ taking care of the 
fact that one sum is over positive index, while the other 
extends over the interval [1 — n,n\. Then, one computes 




■ (fe) ■ (fe) 



for 



for 



J>0 
J <0. 



This linear transformation is shortly noted as: 



(54) 
(55) 



where Ti. is the above detailed transformation. The form 
of Q is dictated by Ea. (|49b ') for positive index, and from 
the fact that the inverse Fourier transform of b (see be- 
low) is real. The new sampling series is finally obtained 
from: 

Jfe+i)^ J7r^^_ij^(,)j_ (56) 

The equations (|50I56|) define one step in the algorithm 
relating u;^''~^^'> to u'^'^K We note this step w^'^^-^^ = 
T{ijj'^''^). The searched function is clearly a fixed point 
of the transformation T defined above in a discretized 
version. The uniqueness of the transformation O results 
from that of the harmonic field in the domain S with an 
equipotential condition on the boundary and a constant 
gradient perpendicular to the boundary at infinite dis- 
tance from it. Therefore, the only condition to consider 
is the stability of the fixed point. As it can be seen in 
Rcf. 1], this fixed point is attractive at the only condi- 
tion that: 



max |/i'(u)| < 1 - 



(57) 



The procedure described above gives a very fast and ef- 
ficient way of building conformal mappings. It has to be 
noted indeed that this two-dimensional problem is solved 
here by the only use of a few one-dimensional fast Fourier 
transforms. In cases of slopes locally overcoming one, it 
is possible to use under-relaxation schemes (see details 
in Ref. 0,0|) or to decompose the mapping in several 
substeps such that the convergence critcrium is always 
fulfilled. 
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